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(54) Method and device for assessing the stability of an electric power transmission network 

(57) In a method and system for assessing the sta- 
bility of an electric power transmission network, where 
at least one pair of measurements comprising a first and 
a second measurement point (P1.P2), each measure- 
ment point comprising a voltage and a current phasor, 
is processed and where a Thevenin impedance (Zt) and 
a present stability margin (dS(k)) value are computed, 
a validity indicator (v) is computed which depends on 
whether there is a difference between the first and sec- 
ond measurement points and whether there is a differ- 
ence between corresponding estimated first and second 
load impedances (Za). From all validity indicators (v) as- 
sociated with all of the at least one pair of measurements 
a quality indicator (q) is computed that is associated with 
the Thevenin impedance value (Zt) and with the present 
stability margin (dS(k)). 

The invention allows to continuously compute and 
output a present stability margin (dS(k)) value and to 
provide a measure (q) of its quality. 
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Background of the Invention 
l0 00 2 , E,e* ri cpowert ra ^^ 

ng geographically separated regions, and ^^^^^SSX^^ * a central energy manage " 
lines. Power generation and load flow in a network wrth several * the network stable, i.e. 

Lntsystem^n important issue in the con,^ 

to avoid voltage collapses and sw.ngs. A method for assessing newer* y ^ Computer 

in the paper "Use of local measurements ^^^^^^^S^^ and its applications", K. Vu et 
Applications (PICA) May 12-16, 1997 IEEE and hn ^^f™^ are herewith incorporated by reference, 
a... Power Systems Computation <^^^^^iJ^«i M . currents and voltages locally in order to 
These articles describe a "Voltage Instab.lrty Mtf (V IP e , Qne part o{ an electric power system 

infer a proximity to voltage col.apse/The <^^™££^\ represented by its Thevenin equrvalent 1 
is treated as a power source, another part as a , oad „ presented by an apparent load imped- 

with a source voltage Et and a ThSyen.n " Za are estimated from the current and 

impedance Za1 is computed as 



1000., aW .„asecon am e as u,.™„,p<>,n,oo m p^n9 a s«»n !l vo« i , 9 .p ha so,V 2a na a seco,dcu ro n,ph aS or 12 , 
the Thevenin impedance Zt is computed as 



v+ V2-V1 

35 Zt= -1^TT 
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[0 005, Note that the phasors as we,, as the impedances are ™"^^ 
L Thevenin impedance Zt correctty .n , prtcuto. n \^^%^Z^ooe Zt is invalid, 
second measurement points does not ^change ,*m the ; a,u ^ 

[0006] In USpatent5,631 ,569, a system for temporary ^ eas " reman ^ ™ n f|uct ; ations in the power network, which 
of the Thevenin impedance Zt is estimated by ^.ng advantage of ^ ^ l^tege ^ an RMS current are meas- 
cause measured vottages and currents to change^ RMS £££ZZZ£^ Tne 9 Th6ve nin impedance is only com- 

falls below one sigma, the mean of the values .s d '^«* occasiona | dis piay but a continuous generation of esti- 
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and a power network stability estimating device according to the claims 1 and 10, respectively, and by a computer 
program product according to claim 11 . 

[001 0] In the inventive method and device for assessing the stability of an electric power transmission network, where 
at least one pair of measurements comprising a first and a second measurement point, each measurement point 

s comprising a voltage and a current phasor, is processed and where a Thevenin impedance and a stability margin value 
are computed, a validity indicator is computed which depends on whether there exists a predetermined difference 
between the first and second measurement points and, if this is the case, also on whether there exists a predetermined 
difference between corresponding first and second load impedances. From all validity indicators associated with all of 
the at least one pair of measurements a quality indicator is computed that is associated with the Thevenin impedance 

10 value and with the stability margin value. 

[0011] The invention allows to continuously compute and output a stability margin value and to provide a measure 
of its quality. This allows an operator to take informed control actions. In a preferred embodiment of the invention, the 
stability margin value and associated quality indicator are transmitted to a network or substation control system. This 
allows, for example, the control system to automatically react to a problem in the network, indicated by the stability 

is margin, only when the associated quality indicator qs show that the stability margin is of a high quality, i.e. higher than 
a predetermined threshold value. 

[0012] In a preferred variant of the invention, the validity indicator v indicates whether the associated pair of meas- 
urements provides a valid, an unchanged or an invalid Thevenin impedance. 

[0013] In a preferred embodiment of the invention, the device comprises means for displaying the quality indicator 
20 q and/or means for transmitting the quality indicator q and the stability margin value to a control system that controls 
at least part of the power network. 

[0014] A computer program product according to the invention comprises a computer readable medium, having 
thereon: computer program code means to make, when said program is loaded, the computer execute a method for 
assessing the stability of an electric power transmission network according to inventive method. 
25 [0015] Further preferred embodiments are evident from the dependent patent claims. 

Brief Description of the Drawings 

[0016] The subject matter of the invention will be explained in more detail in the following text with reference to 
30 preferred exemplary embodiments which are illustrated in the attached drawings, in which: 

Figure 1 shows schematically a conceptual structure for assessing network stability, according to the state of the art; 
Figure 2 is a diagram showing currents and voltages at a selected point in an electric power network; 
Figure 3 is a diagram showing two measurement points in an impedance plane; 
35 Figure 4 is a diagram showing a set of measurement points in the impedance plane; and 
Figure 5 is a flow diagram of the inventive method. 

[0017] The reference symbols used in the drawings, and their meanings, are listed in summary form in the list of 
reference symbols. In principle, identical parts are provided with the same reference symbols in the figures. 



40 



Detailed Description of Preferred Embodiments 



[0018] Figure 1 schematically shows a conceptual structure for assessing network stability, in which one part of an 
electric power system is treated as a power source, another part as a load. The power source is represented by its 

45 Thevenin equivalent 1 with a Thevenin or source impedance Zt. The load is represented by an apparent load impedance 
Za. A VIP device 2 determines phasor data with a phasor measurement unit residing, for example, at a feeder at the 
bay level of substations or at branching points along transmission lines. A voltage phasor represents, for example, a 
voltage of the feeder or line, while a current phasor represents current flowing through the feeder or line. 
[0019] The phasor data represents a phasor and may be a polar number, the absolute value of which corresponds 

50 to either the real magnitude or the RMS value of a quantity, and the phase argument to the phase angle at zero time. 
Alternatively, the phasor may be a complex number having real and imaginary parts or the phasor may use rectangular 
or exponential notation. Phasors may be used to represent quantities such as the voltage, current, power or energy 
associated with a phase conductor or an electronic circuit. By contrast, conventional sensing devices used in power 
networks generally measure only scalar, average representations, such as the RMS value of a voltage, current etc. 

55 [0020] In some VIP applications, the phasor data is collected from phasor measurement units that are distributed 
over a large geographic area, i.e. over tens to hundreds of kilometres. For applications in which the phasor data from 
these disparate sources are analysed in conjunction, they must refer to a common phase reference. In other words, 
the different phasor measurement units must have local clocks that are synchronised to within a given precision. Such 



3 



EP1 211 775 A1 



10 



15 



20 



25 



30 



35 



40 



45 



50 



55 



a synchronisation of the phasor measurement units is preferably achieved with a known time distribution system for 

rSX—ng (GPS) system. In a typical implementation the phasor mH 
everv 200 or every 100 or preferably every 40 milliseconds, with a temporal resolution o preferably less , than 1 mrth 
second Tn a peered embodiment of the invention, the temporal resolution is less than 10 <™croseconds, which cor- 
responds to a ? phase error of 0.2 degrees. Each measurement is associated with a time stamp denved from the syn- 
^hrrvnieoH iftral clock The ohasor data therefore comprises time stamp data. 

raoS Th VP assesses the stability of the electric power transmission network by determining a stab.hty margin 

through a given transmission line of the network. A global power marg.n combines phasor data collected from a plurality 

SSmJbl at a feeder of a tie line or a load. When the load gets stronger, the current increases and the voltage 
drepsThelcll ££££ current curve (V/l curve) is represented by the drawn out line 3. The actual bjhavjour 
fs influenced bv the entire network, but is estimated from local measurements by the VIP dev.ce. The VIP computes 
oroide an esSrTte o this curve, shown by the dashed line 4. In particular, the slope of the curve correspond to the 
magnttude iSSE Thevenin impedance Zt computed by the VIP device. Given this estimated curve the^ maximum 
power in MVA that may be transmitted through said feeder under the current network conditions . > <™p*£ The 
21 power corresponds to a point on the estimated V/l curve 4 for which the product of current and votege « 
mSrCJSSgln. for a present operating point 5, is the difference between 

power being currently transmitted. Let the present time current and voltage magnitudes be l(k and V(k), respectively. 
Then the present power margin or present stability margin dS(k) correspondmg to said feeder is 

10024"' The^cepts of the invention are now explained in terms of a projection into the ^^^ZZZ 
Diane The invention however involves vector values representing phasors for voltages, currents and .mpedances for 
Shich the ^ame concepts hold. Only in the computation of the stability margin dS shown above are scalar values, i.e. 

TooSrSu^ fi-t measurement point P1 and a second measurement point P2 in an 

7SZ Sr C e Ine 7 ShS impedance line 8 which pass through the axes intention and correspond e^ 
toT aZrem load impedance Za ±1 0%. All measurements that differ from the first measurement point P I by less 
than a X maximum vector deference lie in a circle 9 centred at the first measurement point P . In order to compute 
sa^d SencTvottage and current values are scaled, e.g. with respect to nominal operating va ues A source hnped- 
rnr^nS?assing trough both the first measurement point P1 and the second measurement pom :PJ ha a slope 
rhTcorresDonds to the value of the Thevenin impedance Zt computed from the two measurement points P1P2^A 
meas^m^poiS that lies on the first impedance line 6 of the first measurement point P1 does not provide, in con- 
ation with the first measurement point P1 . information about the s.ope^ ^^^^^ m in tne 
K>0261 Fioure 4 is a diagram showing a set of successive measurement points P(k-n) ... P(k-2) , k<k i , kjkj i . 
mpedanc^ plane. Accordig to the invention, any pair of points may be used in an attempt to 

impedance Zt. The invention tests a preferably large number of pairs, determines whK* pairs provide valid values, and 
computes a composite Thevenin impedance Zt along with a measure of its quality. 

mm The inventive method functions as follows: The VIP device 2 repeatedly measures, at a given point ^ ^the 
powe syltem andTa plurality of time instants, a plurality of measurement points, each consisting of a votege phasor 
and a cunL phasor. At least one pair of measurements is determined, comprising a f.rst measurement point P1 and 
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a second measurement point P2, where the first measurement point P1 comprises a first voltage phasor V1 and a first 
current phasor 11 and the second measurement point P2 comprises a second vottage phasor V2 and a second current 
phasor 12. 

[0028] In a first preferred embodiment of the invention, each pair of measurements corresponds to a pair of succes- 
5 sive measurement points. That is, given a sequence of measurement points P(k-n) ... P(k-2) , P(k-1) , P(k), a first pair 
of measurements is{P(k-n), P(k-n+1}, a second pair of measurements is{P(k-n+1), P(k-n+2}, etc... upto{P(k-1), P(k)}. 
[0029] In a second preferred embodiment of the invention, a predetermined number of measurement points is meas- 
ured and stored. A plurality of pairs of measurements is determined by pairing each member of the set of measurement 
points with all other members of the set. That is, given a sequence of n+1 measurement points P(k-n) ... P(k-2), P(k- 
10 1) , P(k), pairs are, with m=k-n: 

{P(m), P(m+1)}, {P(m), P(rn+2)}, {P(m), P(m+3)}, ... {P(m), P(k))} f {P(m+1) ( P(m+2)}, {P(m+1), P(m+3)}, {P(m+1), P 
(m+4)}, ... {P(m+1), P(k)}, {P(m+2), P(m+3)}, {P(m+2), P(m+4)}, {P(m+2), P(m+5)} t ... {P(m+2) ( P(k)}, 

{P(k-3), P(k-2)}, {P(k-3), P(k-1)}, {P(k-3), P(k)}, {P(k-2), P(k-1)} f {P(k-2), P(k)}, {P(k-1), P(k)}. 

is Each ordered pair of measurements, consisting of an older and a newer measurement point, corresponds to the first 
measurement point P1 and second measurement point P2 in the repeated execution of the following steps: For each 
pair of measurements, a first boolean value that indicates a change or difference between the first and second meas- 
urement point is computed. The first boolean value is true only if a minimal difference between the first and second 
measurement point (P1 ,P2) exists. 

20 [0030] In a preferred embodiment of the invention, the first boolean value indicates whether at least one of a voltage 
difference between the first and second voltage phasor and a current difference between the first and second current 
phasor exceed a given respective voltage or current difference threshold. Referring to Figure 3, this corresponds to a 
test whether the second measurement point P2 lies outside a rectangle (not shown) parallel to the axes and centered 
around the first measurement point P1 . In a preferred embodiment of the invention, the voltage threshold is approxi- 

25 mately less than one or one half percent of the operating voltage, and the current threshold is approximately 1 to 2 
percent of the operating current. 

[0031 ] In another preferred embodiment of the invention, the first boolean value indicates whether a vector difference 
between the first and second measurement point exceeds a given distance threshold. Referring to Figure 3, this cor- 
responds to a test whether the second measurement point P2 lies outside the circle 9 or an oval (not shown) centered 
30 around the first measurement point P1 . The exact shape of the circle 9 or oval depends on scaling factors for currents 
and voltages. 

[0032] Since the impedances are vector values, the circle, oval and rectangle referenced in the above explanation 
correspond to their higher-dimensional equivalents. 

[0033] If the first boolean value is false, then there is not a minimal difference between the first and second meas- 
35 urement point P1 ,P2, that is, the second measurement point P2 lies within said rectangle or oval, then the measurement 
points P1 ,P2 are considered to be essentially equal and a corresponding counter in incremented. In the embodiment 
in which each pair of measurements corresponds to a pair of successive measurement points, this equality indicates 
that a network state has not changed essentially since the last valid Thevenin impedance Zt was computed and it may 
therefore be displayed or used by a control system. 
40 [0034] If there is a minimal difference between the first and second measurement point P1 ,P2, then, from the first 
measurement point, a first load impedance 



45 



za,.£. 

[0035] is computed, and, from the second measurement point, a second load impedance 



is computed, and a second boolean value is computed that indicates whether a difference between the first and second 
load impedance exceeds a given impedance threshold. Said difference is preferably expressed in terms of one of the 
load impedances, for example in that the second load impedance should differ from the first load impedance by more 
than ±5%. Referring to Figure 3, this corresponds to the condition that the second measurement point P2 should lie 
55 outside the sector between the second impedance line 7 and the third impedance line 8. 

[0036] If the difference between the first and second load impedance does not exceed the given impedance threshold, 
then the network state has changed, but the pair of measurements cannot be used to provide a valid estimate of the 
Thevenin impedance Zt, and a corresponding counter is incremented. 
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meni and a corresponding counter is incremented. The Thevenin impedance Zt ,s computed as 



^ 12-11 



spending to unchanged, invalid and valid Th6veni 7T^nT G vL 2QQ ms over a duration G f 10 seconds results in 

P *ra»=™S 

S 4 ^ 7 ^^' S'S'r define* how man, pa,„ o, ,ha se. o, pa*, a. 



^w^N^Ws+wS-Nval 
' Ntot 



margin dS(k) at a present time is computed as 



ir ,„ (V(k)-|Zt| ; l(k))2 
dS(k)- ^ 



where V(k) is a magnitude of a present voltage and l(k) is a present current magnftude. both corresponding to the 

^.T^p^^ 

KTS — a m e»,od aad «*. gl va a po—y updated «ab,»» 

^!;.ra:r=t^rr,aS 
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gration of values across the power system. 

[0047] The additional information provided by the quality indicator q allows a differentiated use by either an operator 
or an automated process, as shown in the following examples: 

[0048] If more than half of ail pairs of measurements of a set have essentially equal currents and voltages, then the 
5 power system does not show any dynamic variation and no new values need to be computed. An existing valid value 
of Thevenin impedance Zt or stability margin dS remains valid. 

[0049] If more than half of all pairs of measurements of a set have essentially different currents and voltages, but 
less than half have essentially a different load impedance, then the power system does show a dynamic variation, but 
without large load changes, and a Thevenin impedance Zt or stability margin dS should have only informative character. 
*o [0050] If more than half of all pairs of measurements of a set have essentially different currents and voltages, and 
more than half have essentially a different load impedance, then the power system does show a dynamic variation 
with large load changes. ATh6venin impedance Zt or stability margin dS can be used to decide on protective measures 
such as load shedding. 

[0051 ] In a further preferred embodiment of the invention, depending on the quality indicator and the present stability 

is margin dS(k) value, a control action is initiated automatically. In order to do this, Thevenin impedances Zt and associated 
stability margins dS are transmitted to control system or energy management system of a power network, or to a control 
system of an electric substation. Since certain critical control actions such as load shedding are required if network 
stability is in danger, a corresponding control program effects the shedding of a given load only if the stability margin 
dS has a predetermined minimal quality. For a less important load, a lower quality may be required. 

20 [0052] Figure 5 is a flow diagram of the inventive method. In step 51 , time stamped voltage and current values are 
measured. In step 52, these values are stored in a buffer, providing a set of measurement points on which the subse- 
quent steps operate. In step 53, one pair of measurements is selected and a check if made if at least the voltage or 
the current values of the two measurements differ. If this is not the case, i.e. if both the voltage and current values are 
unchanged, then, in step 54, a counter for the number of measurement pairs with unchanged Thevenin impedance Zt 

25 is incremented. Otherwise, in step 55, it is checked whether he load impedances corresponding to the pair of meas- 
urements differ. If they do not, then, in step 56, a counter for the number of measurement pairs with invalid Thevenin 
impedance Zt is incremented. If they do, then in step 57 the Th6venin impedance Zt is computed and stored, and, in 
step 58, a counter for the number of measurement pairs with valid Thevenin impedance Zt is incremented. After each 
of the incrementing steps 54,56,58, step 59 tests if all measurement pairs have been processed. If this is not the case, 

30 then processing resumes at step 53. Otherwise, in step 60, the quality indicator and best Th6venin impedance estimate 
are computed from the stored Thevenin impedance Zt and from the values of the counters. The counters are reset to 
zero, and in step 61 the best Thevenin impedance estimate and quality indicator are displayed and/or transmitted to 
a control system, together with associated timing information. 

35 List of designations 
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dS / k) present stability margin 

Za apparent load impedance 
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test if all measurement pairs have been processed 
compute quality indicator and best Thevenin impedance estimate 
display and/or transmit results 
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Za1 = -pp. 



g) computing, on.y « the first boo.ean va.ue is true, from the second measurement point (P2) a second .oad 
impedance 



Za2= 1 2 > 
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va'u^ a vaiidtty indicator (v) that is associated with the pa.r of measurements, 
measurements, a composite Thevenin impedance Zt, and 
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k) computing, from all validity indicators (v) that are respectively associated with all of the at least one pair of 
measurements, a quality indicator (q) that is associated with the composite Thevenin impedance value and 
with the present stability margin dS(k). 

2. Method according to claim 1 , characterized in that the stability margin dS(k) at a present time is computed as 

dS(k y (V(kHZt|.l(k))g 

where V(k) is a magnitude of a present voltage and l(k) is a present current magnitude, both corresponding to the 
present measurement point (P(k)). 

3. Method according to claim 1 , characterized in that the first boolean value is true only if at least one of a voltage 
*5 difference between the first and second voltage phasor and a current difference between the first and second 

current phasor exceed a given respective threshold. 

4. Method according to claim 1 , characterized in that the first boolean value is true only if a vector difference between 
the first and second measurement point (P1 ,P2) exceeds a given distance threshold. 

20 

5. Method according to claim 1 , characterized in that the repeated execution of the method steps applies to suc- 
cessive measurement points. 

6. Method according to claim 1 , characterized in that, 

25 in step a), a set of measurement points is measured, and that, 

in step b), a plurality of pairs of measurements is determined by repeatedly pairing each member of the set of 
measurement points with all other members. 

7. Method according to claim 1 , characterized in that the validity indicator (v) indicates whether the associated pair 
30 of measurements provides a valid, an unchanged or an invalid Thevenin impedance. 

8. Method according to claim 1 , characterized in that a visual representation of the quality indicator (q) is displayed 
on a display device. 

35 9. Method according to claim 1 , characterized in that, depending on the quality indicator (q) and the present stability 
margin dS(k), a control action is initiated. 

10. Power network stability estimating device, where said network comprises a plurality of substations, buses and 
lines, where the device is configurable to measuring, at a given point of the power network and at a plurality of 
40 time instants, a plurality of measurement points, each consisting of a voltage phasor and a current phasor, and 

where the device comprises 

a) a means for determining at least one pair of measurements comprising a first measurement point P1 and 
a second measurement point P2, where the first measurement point P1 comprises a first voltage phasor V1 

45 and a first current phasor 11 and the second measurement point P2 comprises a second voltage phasor V2 

and a second current phasor 12, 

b) a means for computing a Thevenin impedance Zt as 

Zt - V 2-V1 

50 a- l^PT' 

c) and a means for computing, from the Thevenin impedance Zt and from a present measurement point (P 
(k)), a present stability margin dS(k) that characterizes network stability, 

characterized in that the device comprises 
55 e) a means for computing, for each of the at least one pair of measurements, a first boolean value that is true 

only if a minimal difference between the first and second measurement point (P1 ,P2) exists, 
f) a means for computing, only if the first boolean value is true, for each of the at least one pair of measurements, 
from the first measurement point, a first load impedance 
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Za1= TP 



g, a means for computing, on.y ii the first boolean va.ue is true, for each of the at least one pair of measure- 
ments, from the second measurement point, a second load .mpedance 



Za2 =T2' 



JiSSL. vi. £«, indoor ,v, M k -~^^£SEEL«-~ 
a composite Thevenin impedance Zt, and f measure . 

20 present stability margin dS(k). 

electric power transmission network according to claim 1 . 
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